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Abstract 

The half-filled extended Hubbard model, in one and two dimensions, is studied by 
means of the 2-pole approximation within the Composite Operator Method with 
the aim at improving the possibilities to describe some of the experimental features 
observed for quasi-lD organic superconductors and Cu-0 planes of cuprates. The 
phase diagrams (T-V and V-U) are analyzed with respect to the paramagnetic metal 
- paramagnetic insulator - charge ordered phase transitions. The relevant features 
of the diagrams (rank of the phase transitions, critical points, reentrant behavior) 
are discussed in detail. 
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1 Introduction 



The extended Hubbard model [1] is the simplest Hamiltonian that takes into 
account non-local Coulomb interactions. Its relevance to the description of 
cuprate superconductors, fullerides and materials presenting charge ordered 
phases (rare-earth pnictides, vanadium oxides, manganites, magnetite, Bech- 
gaard salts, . . . ) is widely recognized in the literature [2,3,4,5]. 

In this paper, we study the model, at half-filling in one and two dimensions, 
by means of the 2-pole approximation within the Composite Operator Method 
(COM) (see Ref. [6] and references therein). The latter, starting from the con- 
sideration that the original electrons lose their identity in presence of strong 
correlations, uses composite operators (i.e., operators built as products of the 
electronic ones) to describe the effective elementary excitations of the inter- 
acting system under analysis. The non-canonical algebra satisfied by the com- 
posite operators dictates constraints (i.e., Algebra Constraints [6]) that are 
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systematically used to fix the representation where the propagators are real- 
ized and self-consistently compute the unknown parameters appearing in the 
theory [6]. 

In the next sections, we will, first, briefly present the model and the method 
and, then, discuss the phase diagrams obtained by the application of the latter 
to the former. The rank of the transitions between homogeneous and charge 
ordered phases are studied together with the relation between such transi- 
tions and the paramagnetic metal-insulator one. The relevant features of the 
diagrams (critical points and reentrant behavior) are also analyzed in detail. 



2 Hamiltonian, field equations and solution 

The extended Hubbard model reads as 



H = Y\-VC ] {i)c{i) - 2dtc\i)c a {i) + Un^n^i) + dVn{i)n a {i)} (1) 



where \i is the chemical potential, c^(i) = (cl(i) c\(i)) is the creation electronic 
operator in spinorial notation, i — (i, t), i is a lattice vector of a rf-dimensional 
square lattice, t is the hopping integral, U is the onsite Coulomb interaction, 
V is the intersite one, n(i) = n^(i) + n±(i), n a (i) is the number operator for 
electrons of spin a. Hereafter, the hopping integral will be used as reference 
unit for all energies. We have used the notation a (i, t) = J2j a y0(j>£)> where 
can be any operator and ay is the projector on the first 2d neighbor sites on 
the lattice. We have a(k) = ^"[ay] = l/dZ)^ =1 cos(A; n ), where T is the Fourier 
transform. 

We have chosen, as basic field, = (i) rf (i)) , where = n(i)c(i) and 
r)(i) = c(i) — = [1 — n(i)]c(i) are the Hubbard operators. They satisfy the 
following equations of motion 



where ir(i) = \a^n^(i)c a (i) + c(i)c a \i)c(i), = (—!,<?), n^i) = (n(i),n(i)) 



is the charge and spin number operator, n(i) = c\i)ac(i) and a are the Pauli 



After the choice we made for the basis, we have computed, in the 2-pole ap- 
proximation within COM [6], the retarded thermal Green's function G(k, oo) = 




-fjg(i) - 2d[tc a (i) + tn(i) - Vi(i)n a (i)} 
K -(H - U)r]{i) + 2d[tir(i) + Vrj{i)n a (i)} } 



(2) 



matrices. 
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J 7 l(R\^{i)il)\j)\) (1Z is the retarded operator) 

where Ei(k) are the eigenvalues of the energy matrix e(k) = m(k)/ _1 (k), 
/(k) = ^F({ip(i,t),il)'(j,t)}) is the normalization matrix of the basis and 
m (k) = J r ({i^(i,t),V t (j> i )})- The spectral weights <jW(k) can be computed 

as 

n 

(7«(k) = X;n«(k)nr 1 (k)/ c6 (k) a,6 = l,...,n (4) 

c=l 

where the matrix fi(k) has the eigenvectors of e(k) as columns [6]. For the 
sake of brevity, we do not report the complete expressions of Ei and 
which can be found in Ref. [7]. The parameters appearing in the theory, 
and not connected to the Green's function under analysis ((n(i)n a (i)) and 

\{ n< u{i) n Lt(i)) ~ (t c t W c iW]° c l W C T W))' wm be self-consistently computed (the 
first) by calculating the density-density correlation function within the one- 
loop approximation [8] and (the second) by means of the local algebra con- 
straint [6] (£{i)r]\i)) = 0. 

The set of self-consistent equations, fixing the parameters appearing in the 
theory, is highly non-linear. According to this, it is not strange at all that it 
admits two distinct solutions that, hereafter, we will call COM1 and COM2. 
The Composite Operator Method tries to give answers in the whole space 
of model and physical parameters and the presence of two solutions should 
be seen as a richness of the method. For the one-dimensional system we will 
consider only COM2 solution as with this, for the simple Hubbard model, we 
got excellent agreements with the Bethe Ansatz exact solution [9]. For the 
two-dimensional case, we will study both solutions as they will permit us to 
analyze two different behaviors that could be both observed experimentally. 



3 Phase diagrams 

3.1 One- dimensional system 

COM2 solution for the one dimensional system reports a phase transition from 
the Mott insulating phase to a inhomogeneous charge ordered state of checker- 
board type for positive values of the intersite Coulomb potential greater than 
some critical one. This result is consistent with many other studies [10,11,12]. 
The nature of this phase transition is not well understood yet and currently 
under intense investigation [11,12]. The phase diagram in the plane V-U is 
shown in Fig. 1 (left panel). The phase transition is of the second order up to 
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Fig. 1. (left) The phase diagram V-U at zero temperature (PI = paramagnetic 
insulator, CO = charge ordered phase), (right) The phase diagram T-V at U = 2, 
4 and 8. 

V ~ 2.5; it becomes of the first order for higher values of V. The kind of phase 
transition that we are here analyzing is completely different from the one usu- 
ally reported in the literature in proximity of the U = IV line as this latter 
occurs between an inhomogeneous charge ordered phase and an homogenous 
spin ordered phase. 

In Fig. 1 (right panel) we report the phase diagram in the plane T-V for 
U — 2, 4 and 8. For U = 2 the transition is second order for all values of T. 
For U = 4 we can see a first signature of a reentrant behavior. By this latter 
we mean a situation in which by increasing the temperature we can first enter 
and then exit a phase when within another. For U = 8 the reentrant behavior 
is clearly evident. It is interesting to observe that the transition is continuous if 
there is no reentrant behavior. When instead there is a reentrant behavior, the 
transition is discontinuous up to the turning point, then becomes continuous. 
It is worth noticing that the transition is clearly marked by a discontinuity in 
the nearest-neighbor density-density correlation function. 



3.2 Two-dimensional system: COM2 



COM2 solution for the two-dimensional case has similar characteristics to that 
observed in the ID case. The phase diagram in the plane T-V is reported in 
Fig. 2. As regards the V-U phase diagram, at zero temperature, we just have 
a transition line with slope ~ 3 and the transition is continuous for U < 1.8 
and first order for higher values of U. For finite temperature and U — 8 a 
reentrant behavior as function of temperature is clearly observed. The transi- 
tion is first order up to the turning point T = 0.6, then becomes continuous. 
For U < 1.8 no reentrant behavior is observed. The fact that charge ordering 
may disappear by decreasing temperature has been experimentally observed in 
Pr . 65 (Ca . 7 Sr . 3 ) . 35 MnO 3 [13] and La 2 _ 2x Sr 1+2x Mn 2 7 (0.47 < x < 0.62) 
[14,15]. 
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Fig. 2. The phase diagram T-V at U = 8. 




Fig. 3. (left) The phase diagram V-U at zero temperature (PM = paramagnetic 
metal), (right) The phase diagram T-V at U = 8 and 15. 

3.3 Two-dimensional system: C0M1 



The system undergoes a metal-insulator transition at some critical value of the 
onsite Coulomb potential U, which depends on the intensity of the intersite 
potential V. With respect to the V = case (the simple Hubbard model), the 
metallic region is compressed by the presence of the intersite interaction and 
disappears for V > 5.7. By further increasing the intersite Coulomb potential, 
the system undergoes a transition to a charge ordered state. The complete 
phase diagram in the plane V-U is shown in Fig. 3 (left panel). The diagram 
is characterized by two critical curves, which separate the different phases. 
The lower one controls the MIT, the upper one controls the transition to a 
charge ordered state. The first transition is first order for U < 12 and second 
order for higher values of U; the second one is first order for U > 1.9 and 
second order for lower value of U . 

The phase diagram in the plane T-V is shown in Fig. 3 (right panel). For 
U — 15 there is no metallic phase and we have a critical temperature where 
a transition from an insulating to a charge ordered state is observed. The 
transition is first order up to T = 0.95, then becomes continuous. A reentrant 
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behavior is observed with the same characteristics previously discussed. For 
U — 8 we have two critical temperatures, which characterize the MIT transi- 
tion and the insulator-charge order transition, respectively. Also in this case 
a reentrant behavior is observed in the latter transition. 

4 Conclusions 

The phase diagrams, in the planes V-U and T-V, of the half-filled extended 
Hubbard model, in one and two dimensions, has been studied by means of the 
2-pole approximation within the Composite Operator Method. Transitions 
between the paramagnetic (metal and insulator) phase and a charge ordered 
state of checkerboard type have been found. Reentrant temperature behavior 
in the plane T-V has been observed with characteristic similar to that exper- 
imentally found for some manganites. The rank of the phase transitions has 
been studied and identified. 
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